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, Abstract. We study £1 resonances, breakup and fusion reactions for weakly bound Be nuclei. The 

^> . absorbing-boundary condition (ABC) is used to describe both the outgoing and incoming boundary 

' conditions. The neutron continuum plays important roles in response and reaction of neutron drip- 

line nuclei. 



> : INTRODUCTION 
oo 

oo 

O I At the drip line, since the separation energy becomes close to zero, theoretical studies 

of nuclear structure and reaction require continuum wave functions. We have recently 
^ \ investigated an efficient and comprehensive method of treating the continuum [? ? ? 

O I ? ? ]. This is practically identical to the one called "Absorbing Boundary Condition 

(ABC) method" in the chemical reaction studies [? ]. The method allows us to calculate 
the continuum wave functions with the outgoing asymptotic behavior of many-body 
O \ systems. 

3 ■ The essential trick for the treatment of the continuum in the ABC method is to allow 

the infinitesimal imaginary part in the Green's function, ie, to be a function of coordinate 
and finite, ie (r) . The £ (r) should be zero in the interacting region and be positive outside 
^ ' the physically relevant region of space. Wave functions obtained using the ABC method 

. are meaningful only in the interacting region. However, this provides all the necessary 

information to obtain the scattering matrix. 

In order to understand how the ABC method is functioning in later applications, let us 
consider the potential scattering of a particle. The scattering wave function is given by 



I V^^^) = |k) + ^_^^./ |k) = |k) + I WctJ). (1) 
The scattering amplitude, f{^), is usually defined by its asymptotic behavior 



but can be also written in a form 

m) = -^^'^'I^IV^^^) = / ^rexp(-/k'-r)V(r)v/W(r), (3) 



where Q. is the direction of k' and |k'| = |k|. Equation © implies that the f{^) can be 

determined by the scattering wave function, i//^^^(r), in the interacting region where 

y 7^ 0. In other words, the \]/^\r) outside the interacting region is not needed to 
determine the scattering properties. This is why we are allowed to add the absorbing 

potential, — ze(r), to the Hamiltonian as long as the the interacting region is 

correctly described. 

In this paper, we present examples of the ABC method. The method well works in the 
time-dependent and time-independent approaches of mean-field and few -body models. 



APPLICATION OF ABSORBING BOUNDARY CONDITION (ABC) 

Time-dependent-Hartree-Fock (TDHF) calculation in the 3D 

coordinate space 

In studies of giant resonances, effects of the continuum has been treated in the 
random-phase approximation (RPA) with Green's function in the coordinate space [? 
] . However, it is very difficult to directly apply the method to deformed nuclei because 
construction of the Green's function becomes a difficult task for the multi-dimensional 
space. We have shown that the ABC method is very useful to treat the electronic 
continuum in deformed systems, such as molecules and clusters [? ]. We have also 
investigated the applicability of the ABC in studies of nuclear response calculations [? 
? ]. In this section, we discuss properties of the giant dipole resonance (GDR) in nuclei 
at the A^^ = Z line, ^Be, and at the neutron drip line, ^'^Be. 

We use the ABC in the time-dependent Hartree-Fock (TDHF) calculations on a three- 
dimensional (3D) coordinate grid. The time evolution of the TDHF state, det{^,(r,f)}, 
is computed by applying the electric dipole (El) field to the Hartree-Fock (HF) ground 
state: 

/^(^),-(r,0= (^-^v2 + yHF(r,O+Vext(r,O-'e(r)^0,(r,O, /=1~A, (4) 

Ven{r,t)=kra!^^{l-T:z)e-'^^d{t), ra=x,y,z, (5) 

where k should be small enough to validate the linear response approximation. We 
calculate the expectation values of the E 1 operator as a function of time, then Fourier 
transforming to get the energy response [???]. Since all frequencies are contained in 
the initial perturbation, the entire energy response can be calculated with a single time 
evolution. 

We use the Skyrme energy functional of EV8 [? ] with the SGII parameter set. For 
the time evolution of the TDHF state, we follow the standard prescription [? ]. The 
model space is a sphere whose radius is 22 fm. The ze(r) is zero in a region of r < 10 
fm, while it is non-zero at r > 10 fm. The TDHF single-particle wave functions are 
discretized on a rectangular mesh in a 3D real space. Now, we can solve Eq. Q together 



with the vanishing boundary condition at r = 22 fm. Time evolution is carried out up to 
r = 30 h/MeV. 

The density distribution of the ground states of ^Be and ^"^Be both possess a prolate 
superdeformed shape. The El oscillator strengths are shown in Fig.[T] Here, we use a 
smoothing parameter of F = 0.5 MeV. The deformation splitting of the GDR peak is as 
large as about 10 MeV. Since the ground-state shape is almost identical between the two 
isotopes, average peak positions are similar. However, the peak width is very different 
between the isotopes. In '^Be at the neutron drip line, the double-peak structure is almost 
smeared out in the total strength (thick line). On the other hand, we observe prominent 
double-peak structure in ^Be. The significant peak broadening in ^^Be may attribute to 
the small neutron separation energy and the Landau damping. 
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FIGURE 1. Calculated E 1 strength in ^'^^Be. Thin solid and dashed lines indicate the response to dipole 
fields parallel and perpendicular to the symmetry axis, respectively. Thick line shows the total strength. 



Calculation of breakup reaction in the coordinate space 

Let us consider a reaction of a projectile, composed of core (C) plus neutron (n), on 
a target nucleus (T). Denoting the projectile-target relative coordinates by R and the 
neutron-core relative coordinates by r, the Hamiltonian of this three-body system with 
the ABC is expressed as 

2 2 

" ^ " ^ V^nc(r) +y„r(r„r) + Vcr(Rcr) - /e(R,r) (6) 

where /i and m are the reduced masses of projectile-target relative motion and neutron- 
core relative motion, respectively. Vnc, VnT, Vct are the interaction potentials of con- 
stituent particles. The ie vanishes in the region where the nuclear interactions are active. 

The wave function may be expressed as a sum of the Coulomb wave, (R), in the 
incident channel and the scattered wave. 

^(+)(R,r) = v/(+)(R)0o(r) +^scat(R,r) (7) 



where ^o{r) is the ground state of the projectile, described as a n-C bound state. The 
^scat satisfies the following inhomogeneous equation, 

{£ + eo-//}^scat(R,r) = {Kr(r„r) +Vcr(Rcr) -Vc(R)} V/(+)(R)(^o(r), (8) 

where Vc is the Coulomb distorting potential, E is the bombarding energy, and eo is the 
ground-state energy of the projectile. One should note that the right hand side of Eq. ^ 
is a localized function if we can neglect the difference between Vct and Vc at large R. In 
[? ], we have studied a deuteron breakup reaction and compared our results with those 
of the continuum discretized coupled channel (CDCC) calculation. Readers may refer 
to [? ] for numerical details. In this paper, we report the application to a nuclear breakup 
reaction of ^^Be on ^^C. 

The ^^Be-n potential is taken as a Woods-Saxon shape whose depth is set so as to 
produce the 2s orbital binding energy. We adopt the optical potentials for '*^Be-^^C and 
for n-^^C. The radial region up to 30 fm and 50 fm are used for R and r, respectively. 
The /e(i?, r) is non-zero in the region of 20 fm < 7? < 30 fm and in that of 25 fm < r < 

50 fm. The n-^^Be relative angular momenta are included up to / = 3. 

In the left panel of Fig. El we show the elastic breakup cross sections of ^^Be-^^C 
reaction. The filled circles are the result of the ABC calculation and the open circles 
for the eikonal calculation. The elastic breakup cross section is substantially larger than 
that in the eikonal approximation at lower incident energy. The failure of the eikonal 
approximation is apparent at the incident energy below 50 MeV/A. There, the quantum- 
mechanical treatment is required for the three-body continuum. 
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FIGURE 2. (Left) The elastic breakup cross section in "Be-'^C reaction. Quantum calculation with 
ABC (closed circles) is compared with the eikonal calculation (open). 

(Right) Calculated fusion probability for ^ ' Be+^^^Pb as a function of the projectile-target incident energy 
for head-on case (J = 0). Different truncation for n-C partial waves is taken for a dotted line (I < 2) and 
for a dashed line (I < 70). The solid line is the fusion probability for '"Be+^^^^Pb without a halo neutron. 



Time-dependent- wave-packet calculation for fusion probability 



So far, we have used the ABC to simulate the outgoing boundary condition, ie being 
non-zero outside the interacting region. In this section, we set the ie non-zero inside the 



Coulomb barrier, to simulate the incoming boundary condition. 

We consider a fusion reaction of '^Be on ^ospj, j^g^j. ^j^g Coulomb barrier energy. 
Again, ^^Be is described as a weakly-bound system of a neutron and ^^Be core. The 
target is now ^*^^Pb. All the nuclear potentials are taken to be real and of the Woods- 
Saxon type. The ie{RcT) is also of the Woods-Saxon shape and corresponds to the 
imaginary part of the C-T optical potential which describes the fusion between core 
and target. See [? ] for details. The time-dependent Schrodinger equation, 

/^mr,0=//mr,0, (9) 

where H is given by Eq. Q, is solved with an initial wave function 

^(R,r,? = 0) = exp(-/A'i?)$o(R)0o(r), (10) 

where 4>o(R) is a Gaussian wave packet. The initial wave packet, ^»o(R)0o(r), exists 
in a region beyond the range of nuclear interaction of VnT and Vct- The wave packet is 
going to approach to the target, then if the ^*^Be core penetrates the Coulomb barrier, 
the wave function disappears because of absorption by the imaginary potential ie{RcT)- 
The ratio of the flux loss against the initial flux gives the fusion probability. In this 
calculation, since the final destination of neutron is irrelevant, the calculated fusion 
probability contains both complete and incomplete fusions. 

The energy projection [? ] leads to the fusion probability shown in Fig. |21 (right). 
When we include the n-C partial angular momenta only up to / = 2, the result indicates 
a strong fusion enhancement at sub-barrier energies, compared to the two-body fusion 
calculation of ^^Be and ^'^^Pb. However, when we enlarge our model space to include 
/ < 70, this enhancement disappears and we even have suppression. This is consistent 
with similar study for fusion suppression obtained for ^^Be and ^''Ca [? ]. The results 
suggest necessity of high partial waves for the n-C relative motion to describe correct 
dynamics of neutron breakup and transfer under the strong Coulomb field. 



